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Abstract 

According to the recent experimental data of a GSI-experiment, the rate of the number of daughter 
ions 140 Ce 58+ , produced by a nuclear K-shell electron capture (EC) of the H-like ion 140 p r 58 + ) j s 
modulated in time with a period Td = (7.06 ± 8) seconds. We explain this phenomenon by neutrino- 
flavour mixing and show that this can be understood within standard quantum field theory and derive 
a value for the squared mass difference Am^ = m| — m\ = (0.763 ±8) ■ 10~ 4 eV 2 . This proves that 
such processes provide a precise method to investigate neutrino-flavour mixing. 

PACS: 12.15.Ff, 13.15.+g, 23.40.Bw, 26.65,+t 



1 Introduction 



Recently, K-electron capture (EC) and (5 + de- 
cays of H-like 140 Pr 58 + ions were measured in 
the Experimental Storage Ring (ESR) at GSI in 
Darmstadt pQ. It is very exciting that in this 
electron capture process periodic modulations of 
the expected exponential decrease of the number 
of EC-decays per time was observed [2J. The 
rate dNf c (t)/dt of the number of daughter ions 
140 Ce 58+ produced in EC-decay 140 p r 58+ -^ 140 
Ce 58+ + v shows a modulation in time with a pe- 
riod T d = (7.06 ±8) seconds [2J. This value of T d 



is the result of a fit with the function 



dNf c (t) 



dt 



Aec(*) e~ M N m (0) 



(1) 



where iV m (0) is the number of mother ions at the 
production time t = 0, the decay constant A in- 
cludes loss of particles and (3 + decay. The time 
dependent "decay factor" Aec (t) , one cannot call 
it anymore decay constant, was fitted by 

Aec (t) = A EC [1 + a • cos(^ + 0)] (2) 

with a result of the fit, = 27r/T d = 0.89(1) s" 1 
and a = 0.15(3). 
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It was shown in [3j that the K-electron capture 
to positron emission ratios of H-like 140 p r 58 + and 
the He-like 140 p r 57+ ions can be described within 
the standard theory of weak interactions of heavy 
ions [4] with an accuracy better than 3%. 

In recent years neutrino oscillation predicted 
by Bruno Pontccorvo [5] have witness great 
progress both experimentally and theoretically 
[5J, [7]. During the last decades experimental ev- 
idence for neutrino oscillations came from solar 
neutrino experiments (Homestake, Kamiokandc, 
SAGE, GALLEX-GNO, SuperKamiokande (SK) 
and SNO) long-baseline reactor neutrino exper- 
iment (KamLAND), atmospheric neutrino ex- 
periments (SK, MACRO, and Soudan-2) and 
long-baseline accelerator neutrino experiment 
(K2K)i. 

In this article we would like to support the idea 
that neutrino oscillations can also be detected in 
electron capture processes. We show that the ob- 
served time-modulation of the number of daugh- 
ter nuclei in the EC-decay of the H-like 14 °p r 58+ 
ions can be explained by neutrino-flavour mixing 
[8j. Therefore, EC-decays provide a new method 
for studying neutrino mixing and mass differences 
of neutrino flavours. 

In the following we will first discuss the kine- 
matics of the process, then we will start the quan- 
tum theoretical calculations with the calculation 
of the transition matrix clement and finally we 
will determine the decay rate by the methods of 
quantum field theory These calculations should 
prove that we have picked up the right kinemat- 
ics. In the main part of the text the most impor- 
tant steps of the calculations are shownQ After 
these investigations we discuss the relations for 
the amplitude of the oscillations as well as some 
features of the results. 



2 Kinematics 

In the K-capture the mother ion 140 Pr 58+ makes a 
transition to the daughter nucleus 140 Ce 58+ with 

1 The detailed quantum theoretical calculations are shifted 



the mass M d = 130319.252 MeV and the Q-value 

Q = M m - M d = 3348 ± 6 keV. (3) 

In the rest frame of the mother ion 140 p r 58 + W c 
are allowed to take the non-relativistic expres- 
sions for the energy of the daughter nucleus. If a 
neutrino species of mass rrij is emitted we expect 
due to energy and momentum conservation the 
relation 

j^+S- d = Q > (4) 

where kj is the momentum of the neutrino in the 
mass cigenstate rrij. Inserting into this equation 
j = 1 and 2, we expect from the best-fit 

Am| n = ml - m\ = 8.0t° j • 1CT 5 cV 2 (5) 

of the global analysis of the solar-neutrino and 
KamLAND experimental data [6J (see also [9]) 
momenta of k\ w k 2 ~ Q = 3.388 • 10 6 eV and a 
tiny momentum difference of 

ki -ka = 1.12- Kr n cV, (6) 

which can be estimated to be about k\ — k 2 ~ 
Amf 2 /(2<3). For the further kinematic calcula- 
tion we use two neutrino species only. The quan- 
tum theoretical part we will do with a general 
set of neutrino mass eigenstates. In the stan- 
dard procedure for the determination of the de- 
cay rate one integrates over the momenta q of the 
daughter nuclei with mass M d or the momenta 
kj of the emitted neutrino with mass mj . After 
some approximations this could lead to a guess 
of the oscillation period T caused by the mixing 
of the neutrino mass eigenstates m\ and m 2 of 
2i 2 = 27r/(fci - k 2 ) ~ 2tt • 2Q/Am 2 2 = 0.3ms 
which is far away from the experimental value, 
T d ~ 7 s and is disproved by the experiment. 

One can get much closer to the experimen- 
tal value of T d using an oscillation frequency 
ill 2 = Am \ 2 /{2M d ) which one can guess charac- 
terising the emitted neutrinos by their energy E u 
and using energy and momentum conservation for 
neutrinos of different mass and the same energy 
E v . This idea leads to the energy and momentum 

to the appendix. 
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conservation relation E v + (El — m^)/(2M d ) = Q 
and to the prediction T i2 = 13 s [2j- But this idea 
suffers from the problem how to integrate over the 
final neutrino states. In the standard procedure 
one has to integrate over the momenta of the fi- 
nal states, but the two neutrinos with the same 
energy and with different masses which should 
interfere have two phase spaces of different size. 
Which size does one have to use for the integra- 
tion? Even if the momentum difference © is tiny 
compared to the absolute values ki ~ &2 ~ Q of 
the momenta, this problem has to be solved. 
This difficulty with the definition of the phase 



space shows that two involved mass-eigenstates, 
a lighter rrij and a heavier m;, need exactly the 
same momentum k in the creation process when 
they are emitted from the mother nucleus. At 
the same time the emitted neutrinos have to have 
those momenta kj and ki and energies which fol- 
low from the energy and momentum conservation 
relations ^ in the centre of mass system of the 
mother ion. In classical mechanics this looks like 
a contradiction. In quantum mechanics this is 
possible, since the invention of Schrodingers cats 
we are used to have in quantum mechanics ob- 
jects in a superposition of different states. 




Figure 1: Scheme of the momenta involved in the kinematics. The momentum uncertainty k of the mother 
nucleus and the momentum fc; of the neutrino mass eigenstate with mass mi add up to the momentum 
kj of the neutrino with mass rrij. The order of magnitude of the radius kj of the sphere is by many 
magnitudes larger than the magnitude of — k. Therefore the sphere with the radius kj is rather a plane. 
A is the cross-section of this plane with the shaded region of allowed k values. If A shrinks to zero there 
are no oscillations. 



For the two neutrino mass eigenstates the mo- 
menta in the centre of mass system differ by a 
momentum k = kj — ki, sec Fig. [TJ According 
to Eq. (UJ) the lighter neutrino species has the 
larger momentum. Due to the experimental ar- 
rangement in the GSI experiment [2] with ions 
in a storage ring the mother ions have some mo- 
mentum uncertainty and we have to assume a 



wave-packet for the mother ions with a width 
a K for the distribution of the total momentum 
—K. Then neutrinos can be emitted from any 
component of this wave-packet moving against 
the centre of mass system So. 

For simplicity we assume that the mass eigen- 
state rrij is emitted in Sq with momentum k = kj 
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conserving energy and momentum (j4]). To get 
the same momentum for m; in the emission pro- 
cess we have to emit it from a component of 
the mother nucleus which moves with momen- 
tum — k — —(kj — fc;) backwards. In this system 
S- K energy and momentum are conserved for the 
emission of mi (at the vertex of the Feynman di- 
agram) . 

Then, in the centre of mass system So the en- 
ergy conservation is violated for the emission of 
mi, though mi is moving with momentum k\. It 
is the daughter nucleus which moves too fast, 
with momentum — (fc; + 2k). The energy conser- 
vation is violated by a small amount AEi. The 
amount of this violation gives a period Tji of the 
oscillation of the number of decaying nuclei by 
T fl = 2ir/AEi. 

The value of the energy violation AEi for the 
neutrino with mass mi as a function of the neu- 
trino momentum k, the momentum k of the 
mother ion and the angle (p, see Fig. 1, between 
the two momenta is given by 



AEi(k, K,(p) 



(7) 



= W (fc - k) 2 + mf + 



(fc 



2M d 



-Q- 



2M„ 



With Eqs. (@| and ([7]) we can predict from the os- 
cillation period Td in the GSI-experiment the dif- 
ference of the squared masses of the correspond- 
ing neutrinos. Due to the very different scales 
involved in this process the calculations have to 
be done with an accuracy of at least 25 digits. 
The scales range from the mass Md in the order 
of 130 GeV to the energy violation, in the order 
of 6 • 10 _16 cV, i.e. over 27 orders of magnitude. 

For simplicity we start the calculation with ip 
I i and adjust the oscillation period Td — 7.06(8) s 
of the EC-decay of the H-like ion 140 p r 58+ . We 
get as a prediction of the squared mass differ- 



ence Amlj 



7.63(8) • 10~ 5 eV^ 



This value is in good agreement with the value 
• 10~ 5 eV 2 , obtained as a best-fit 



Am\ x 



8.0 



+0.6 
-0.5 



of the global analysis of the solar-neutrino and 
KamLAND experimental data [6] (see also [9]). 
This agreement shows that we are on the right 



track with our picture of the kinematics. 

It still remains to check whether the prediction 
depends on the angle <p, see Fig. 1, between k and 
ki. It turns out that for given Am\ x the value of 
the energy violation AE is very stable against 
variations of tp until close to ±7r/2 with an ac- 
curacy of at least five digits. For ip values larger 
than ir/2 the solution starts to escape. This be- 
haviour is understandable from the geometry and 
the kinematics. Due to the tiny difference k of kj 
and ki the sphere with the radius kj is more a 
plain than a sphere. If <p gets larger than tt/2 
the value of k increases abruptly and within a 
few thousandth of 1° there is no solution of the 
equations possible. Furthermore, k enters mainly 
as k 1 1 , as the component of k parallel to fc; . The 

component normal to ki enters only by the 
term K 1 L /{2M d ) which is negligible small. 

Further, K 2 L /(2M m ) is a negligible small en- 
ergy. Therefore the energy violation is almost 
completely ascribable to the daughter nucleus and 
can be approximated with good accuracy by 



AEi(k,K,p) = 



(fc + k) 2 - k 2 
2M d 



(8) 



where kj and fc; are solutions of Eq. (0| for 
the neutrino mass eigenstates mj and mi and 
k = kj — fc; . For given Am^ the absolute values 
of mj and m; have a negligible influence. There 
is some dependence on the Q-value, but this is 
relatively weak. 

3 Quantum field theoretical 
calculations 

After the discussion of the mechanism we can 
start to use standard quantum field theory to ver- 
ify the applicability of this mechanism. According 
to the standard theory of weak interactions one 
can describe the electron capture by the Hamil- 
tonian 

GpVud 



H w (t) 



V2 



d 6 x [i) n {x)^{l - g A lb)^ P {x)] 



X [$vA x )7li{ 1 ~ lb)i>e{x)], 



(9) 
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where x = (t,r), Gf is the Fermi constant, V u d 
is a CKM matrix element, gA is the axial cou- 
pling constant [6], ijj n (x), ip P {x), ipe{x) and ij) Ut ,(x) 
are operators of the neutron, proton, electron and 
neutrino fields, respectively. 

We take into account neutrino flavour mixing 

m 

i«W= E uti^M ( 10 ) 

j=l,2,3 

by the unitary neutrino-flavour mixing matrix U 

with U e i = C0ST?12 COS $13, t/ e 2 = SUl 1?i2 COS $13 

and U e 3 = sin $13 e~ lS ° F , where are the mix- 
ing angles of the neutrino species j and I. The 
phase Sep is a CP-violating phase [5] which can 
be chosen as e l5cp = ±1. In our analysis the neu- 
trinos Vj (j = 1, 2, 3) are Dirac particles with the 
masses irij and e lScp = 1. 

First we evaluate the matrix element 

(f\H w (t)\i MF ) (11) 

between the initial state (fT3|) of the H-like ion 
i4o Pr 58+ in the state with / = i ( p = i an d 

Mp = ±i and vanishing centre of mass momen- 
tum and the final state 

| /) = |140 Ce 58 +(g> . ( £ ))t) (12) 

of the nucleus 140 Ce 58+ in a + -state with mo- 
mentum q and neutrino mass eigenstate Uj with 
momentum k at time i. According to the dis- 
cussion above we have to allow for the state of 
the mother ions has a more complicated structure 
then for the states of the daughter nuclei. We as- 
sume for the initial state \iu F ) of the mother ions 
in the centre of mass system a wave-packet, cen- 
tred with a width a K around momentum zero 

i ■ \ f d^K ^ — 2" 

l * M * } = J (4^)3/4° " (13) 
x | 140 Pr 58+ («) MF ,i), 

where Mp is the magnetic quantum number of 
the mother ion, with a total angular momentum 
I m = 1. A finite value of a K implies that the 



wave function of the mother nucleus is restricted 
to some region of the centre of mass coordinate 
R. The wave function (|13p is normalised to one 
mother ion. This ansatz is justified by the energy 
uncertainty, introduced by the time differential 
detection of the mother and daughter ions with a 
time resolution shorter than the modulation pe- 
riod [2]. In the final states it is not necessary 
to introduce wave-packets [TU]. For the deter- 
mination of the decay rate one has to integrate 
over all final states and there should be no differ- 
ence whether the final states are counted as plane 
waves or as wave-packets. 

It is useful to treat centre of mass motion, rel- 
ative motion, spin wave functions and time de- 
pendence in the matrix element (f\Hw{t)\iM F ) 
separately In the centre of mass system of the 
wave-packet of the mother ion the wave func- 
tions depending on the centre of mass coordinate 
R read 

f d 3 K -j^T \(R-q-k)R / 14 x 

J V (4^2)3/4 ' U4j 

where the wave functions of neutrino and daugh- 
ter nucleus are normalised to one particle in the 
volume V. Since the wave-packet of the mother 
ion is already normalised (|13p , no further normal- 
isation must be included. 

Due to the point-like nature of the Fermi in- 
teraction the transformation of the electron in a 
neutrino and of a proton in a neutron must take 
place at the same relative coordinate f and re- 
sults in the overlap of the spatial wave functions 
of electron, neutrino and the matrix element of 
the "hadronic" isospin operator t_ between ini- 
tial and final nucleus 

( lia Cc 5S+ (q) Vj(k)\t h _t l + \ lia Pr 5S+ (d)) = 

= : M EC . (15) 

The transition of the electron to the electron neu- 
trino is taken into account with the weak (lep- 
tonic) isospin operator t\. Since the K-electron 
capture is an allowed transition the exponentials 
for the relative wave functions are substituted by 
unity and the matrix clement Mec does not de- 
pend on the momenta k and q. 
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The spinorial part of the matrix element 
(f\Hw{t)\iM F ) we evaluate in the non-relativistic 
limit for neutrons, protons, electrons and mass- 
less neutrino states, with the Dirac spinors 



<P- 



y/2 \-<P 



J' 




(16) 



Here we have omitted additional normalisations 
which would finally cancel with the normalisa- 
tions of the creation and annihilation operators 
and the phase space integrals and assume that 
the spin wave functions ip are normalised to unity. 
Inserting these spinors in the Hamiltonian @ one 
gets within a few lines for the spinor part the well- 
known Fermi- and Gamov- Teller contributions 



[u„(x)7 M (l - gAl^)u p {x)][u u: j{x)-i ll (l - 7 5 )w e (x)] 
= s/2[((p n cp-\tp p tp e ) - gAiVntf-l^^lVpPe)], 

(17) 



where a > acts in the hadronic and leptonic spin 
space, respectively. 

The spin states I m and Id of mother and daugh- 
ter nuclei we decompose in the spin I c of a core 
and proton and neutron spins j Pt „ and include 
angular momentum conservation by 




F = L 



(18) 



where s e and s v are the spins of electron and neu- 
trino. 

To simplify the calculations we quantise the 
spins in the direction of the neutrino momentum 
—q. Therefore, the magnetic quantum number of 
the neutrino state (p~6|) is \i v = —1/2 and due to 
the vanishing spin of the daughter nuclei Id = 



we get Mp = —1/2 and I c = 1/2. In the initial 
state and in the final state we have to couple three 
spin-l/2-states to F = 1 and Mp = —1/2. Since 
the mother nucleus has I m = 1, the initial spin 
state \4>i) is a mixed symmetric doublet state, it 
reads in the representation with magnetic quan- 
tum numbers \fj, c fj, p fj, e ) 



\d>, 



1 



V6 



I tu + m -2 lit). 



(19) 



Due to Id = the final spin state \4>f) is a mixed 
antisymmetric doublet state, it reads in the rep- 
resentation IfJ-cfJ-n^v) 



1 



\4>f) = -7=1 tu - m>. 



(20) 



Since \<f>i) and \<t>t) arc orthogonal states, the 
Fermi contribution in (|17p vanishes. Inserting 
these states in the GT-contribution in p?|) the 
spin operators a ' leave the spin wave function 
of the core untouched, and the spinor part of the 
matrix element (f\Hw(t)\tM F ) reads after a short 
calculation 

- V2gA((pnP-\a h a l \ip p ip e )} =-V&gA- (21) 

We collect now all the terms contributing to the 
matrix element (f\Hw(t)\iM F ) , the spinor con- 
tribution (j2Tj) . the prefactor GpV u d/V% in the 
Hamiltonian ©, the integral over the spatial 
wave functions p5]) , the integral over the centre of 
mass coordinate (|14p , the neutrino mixing pop in 
the wave function, the exponential for the time- 
dependence according to Eq. ([7]) and get with the 
abbreviation 

M GT = -VSg A GpV ud M EC (22) 
the matrix clement 



(f\H w (t)\i MF 



iM, 



GT 



d 3 K —St 



l/(47T 3 a2)3/4 



d 3 R eK*-?-*)* 



U ej ex.p{iAEjt}. (23) 
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For the further calculation it is important to eval- 
uate the R integral carefully for finite extents L, 
or Li, of the volume V and not just to substitute 
the integral by a delta-function. In Eq. (f2"3")) we 
have also taken into account that for Mp = 1/2 
there is no transition possible. 

Since the interaction is weak it is sufficient to 
use time dependent perturbation theory of first 
order to determine the decay rate for the electron 
capture reaction. We average over the two ini- 
tial states \im f ) with Mp = ±1/2, integrate the 
square of the time-integrated matrix element (1231) 
over the possible final states with the density 
(2ir) 3 /V in momentum space and get for the de- 
cay factor 



d 3 fc 



d 3 q 



(2ir) 3 /V J (2ir) 3 /V 



dr(f\H w (T)\i Mf 



(24) 



Here we have inserted an important modification 
of the usual expression which will turn out to be 
necessary for the further calculation. The decay 
rate is derived from the percentage of decayed nu- 
clei by a time derivative and not by a division by 
the time t, as usual. 



Inserting the matrix element (|23|) we can easily 
simplify the expression (j24|) for the decay factor. 
We sum over Mp , integrate the coordinate R over 
the volume V = L 3 and cancel the volume factors. 

For the evaluation of the cf-integral over the 
product of strongly oscillating sin-functions in 
the decay factor (12411 we will use a formula, si- 
miliar to the well-known limit 



dk 



sin 2 (kt) 



f{k) -> Tit / dkf(k)5(k) (25) 



which is valid for large t and for smooth functions 
f(k). The limit 



^sin[(fc - ki)t] sin[(fc - k 2 )t] \k 1 -k 2 \t-»2n 



k — k\ 



k - k 2 



- / dk 



{S(k-h) 



sin[(fc — k 2 )t] 



k - k 2 

sin[(fc — ki)t] 



(26) 



k-k 



k -iB S (k-k 2 )\f(k). 
i > 



is valid for — ^2|i ^> 2n and is discussed in 
the appendix. For a K 3> 1/L we can assume that 
the Gaussians are smooth functions. With these 
calculations we get 



AecW 



d 3 t 



d 3 k 



(2tt) 3 7 (47r 3 a2)3/ 2 



n[ 



2 i II 



i[J2 U ^ [ drejcp{iA^(A!,« J ^}][5^C^ f drexp{-iA^(fc, K ',^)} 



(27) 



The time integrals in Eq. (j2~7| integrate over 
the off-shell contributions J7]| 

ft sin AEjt ae, 

/ drexp{iA^T}= Ag 2 e'— *. (28) 

The sum over the flavour contributions j in 
Eq. (|2"T)) has to be squared and contains inter- 
ference terms. These terms oscillate in time with 



constant amplitude. A division by t, as it is usu- 
ally applied in the determination of the decay rate 
from the square of the amplitude, would destroy 
them. We use therefore the correct definition of 
the rate as the derivative of the square of the am- 
plitude and get for real flavour mixing coefficients 
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u, 



sin— J- ■ Ag Jf i 2 

. - £ p 1 2 

dt\^ e3 *Zl 

j 2 



(29) 



E/£ + J! U ej U e ico S (AEit) 



with interference terms which survive for large t. 
The energy conserving delta-function S(AEj) has 
a negligible dependence on rrij 



dAEj(k,K,ip) 
9m, 



(30) 



and can be moved out of the sum over the neu- 
trino species, ^2jS(AEj) — > S(AE)J2j. Due to 
the smallness of the recoil of the daughter nu- 
cleus the derivatives of AE(k, k, ip) with respect 
to k and k are unmcasurable close to unity. From 
Eq. we get 

dAEAk, k. tp) Q 

J \' '-f-t = 1 + — = 1.000026, 

ok Md 

dAEj(k,K,ip) Q , . 

J \ ' = l--gL = 0.999974. 31 

dn M d v ' 

Therefore, S(AE) selects the modulus k of the 
neutrino momentum in the expression (|27[) for the 
decay factor and modifies its value by some ppm 
only. 

For the further calculations it is simpler to treat 
diagonal and off-diagonal terms separately. In 
the diagonal terms of Eq. (|27|) one can integrate 



the sin function over k[ enforcing k 2 ; = k,[. With 
the help of Eq. (|29|) and using the unitarity of 
the flavour mixing matrix we can drastically sim- 
plify the contribution of the diagonal terms to the 
decay amplitude 



A 



EC 



2 

(7^)3/2 



d 3 fc 
{2-Kf 



(32) 



'6(AE(k,K,ip)). 



Small variations of the centre of mass momentum 
lead to a tiny variation of the momenta of the 
emitted neutrinos. Their contributions to the av- 
eraged decay factor Aec can be summed up by an 
integration over the momentum k of the wave- 
packet of the mother nucleus. The energy con- 
servation S(AE) restricts the momenta k of the 
neutrinos to k = Q. Therefore we get the average 
electron capture constant 



A 



EC 



1« 



2tt 



2- 



[9AG F V ud M EC Qf 



(33) 



The terms in Eq. (|27[) . which are non-diagonal 
in the neutrino flavours, describe the interference 
of the mass eigenstates rrij and mi of the neutri- 
nos with momentum k and lead to an interesting 
oscillation of the decay factor with the frequency 
ft given by the energy violation il = AE. With 
Eq. such a term reads 



A E c« 



2 

x U e 



M GT 



d 3 k yr \f dtti f dn'i 



2a'i 



sm(- 



(2t)3llL(4,3 ff 2)l/2 



Uel 7u f drexp{iAEj(k,K',(p')t} [ dr exp{-iAE t t(k, k, tp)} 
d * Jo Jo 



(34) 



d 3 fc 
(2tt) 3 



n[ 



dKi sin(^-L) 



UejUetfwSiAEjik, 0, tp)) cos(AEi(k, K, <p)t) 
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Here we have used that the integrand is a func- 
tion of k — k' only and that folding two nor- 
malised Gaussian distributions of width a gives 
a normalised Gaussian of width cry/2. In Eq. (|M|) 
both neutrinos, rrij and to/, are integrated over 
the same momenta k. The neutrino rrij is on- 
shell kj = k and the emission of the neutrino m/ 
due to the neutrino momentum ki = k — k in 
the centre of mass system is slightly off-shell by 
AEi(k,K,ip) according to Eq. (0. The value of 
the energy violation AEi(k, K,ip) defines the os- 
cillation frequency in Eq. It's the energy 
violating contributions to the decay rate which 
lead to measurable consequences, to the modula- 
tion of the decay rate in time. 



in Eq. (|34|) leads to a delta- function <5(kj). The 
vanishing momentum k does not allow to solve 
the kinematic equation as can also be seen in 
Fig. 1 where the shaded sphere of allowed ft values 
shrinks to zero and the area A of the overlap with 
the sphere defined by the momentum k = kj gets 
also zero. As a consequence there are no oscilla- 
tions. There are also no oscillations if a K is very 
small. Only if the allowed momenta K are larger 



than \kj — ki \ the sphere with radius k 



kj crosses 



the « cloud in a finite area A which we can esti- 
mate with the minimum of cr\ — (k\ — k<if and 
(27T/L) 2 - (h - fc 2 ) 2 , 

A(\h-k 2 \,L,a K )= (35) 
= Mm{<7 2 - {kt - fc 2 ) 2 , (2tt/L) 2 - (fe x - k 2 ) 2 } 



4 Discussion and concluding 
remarks 

There are two restrictions for the region of al- 
lowed k values, the width er K of the wave-packet 
and the spatial extent L of the wave function of 
the mother ion. If the form of the wave-packet 
is induced by the spatial extent L of the sys- 
tem then there should be a relation of the form 
cr K = 2n/L. In the limit L — > oo the sin- function 



Here we have assumed for simplicity spherical 
symmetry. For an estimate of the oscillation am- 
plitude we have to multiply A with the integrals 
over Ki in the square brackets on the l.h.s. of 
Eq. (|34|) . The product of these three integrals is 
a function of F^{L ■ a K ). This function 



F n {L- 



n 

n[ 

i=l 



is plotted in Fig. 2. 



dfti sin(^L) 

27 ¥ 



icrl 



(36) 




Figure 2: Function F 2 (L ■ cr K ) (left) and function F^(L ■ cr K ) (right), as defined in Eq. (13611 
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Ag C (i) differs from Aec by the factors 
AF^UejUei cos(AEit). Under the assumption 
that the two lightest neutrino flavours contribute 
only we have U e \U e2 = 1/2 ■ sin(2z?i 2 ). We get 
therefore a relation for the oscillation amplitude, 

a = M\h-k 2 \,L^)F 3 (L-a K ) 

(37) 

For simplicity we have assumed that the vol- 
ume V has the same extent in 3 directions. But 
most likely due to the beam geometry a cylindri- 
cal symmetry is more likely with L x sa L y and 
L z = oo. This would mean that the neutrinos 
can interfere only if they are emitted in a direc- 
tion perpendicular to the beam. The product of 
integrals should then be extended over the x and 
y directions only and the function F 2 should be 
used instead of F 3 . 

For the experimental value of the flavour mix- 
ing angle sin(2tf 12 ) = 0.925(32) [6], [9] and the ex- 
perimental value[2] of the amplitude a = 0.15(3) 
we get from Eq. (37]) for the factor A ■ F = 0.32. 

There remains the question why the value 
for the oscillation frequency is so close to the 
value expected from the experimental value of 
the squared mass difference Am^ = m\ — m\ = 
(0.763 ± 8) • 10~ 4 cV 2 j6]. It seems obvious that 
a Lorentz factor for the transformation from the 
centre of mass system to the laboratory sys- 
tem would be necessary. This factor should be 
7 = 1.43 and would increase the oscillation pe- 
riod. Since no essential approximation was made 
in the quantum field theoretical determination of 
this period one has to think about the possible 
reasons for this fact. It seems that we have to 
take into account that the ions don't move in an 
inertial system but around a circle and that the 
emission time of the neutrino takes longer then 
the rotation period. How could the neutrinos lead 
to oscillations of such a long period if the emis- 
sion process would be shorter than the revolution 
period. 

Finally, we would like to give a recipe to pre- 
dict oscillation frequencies in other EC-capture 
systems. First one solves Eq. ([JJ to get k\ and k 2 . 



The mass rrij of the lighter neutrino is not impor- 
tant, as long as it is in a reasonable region, since 
k\ — k 2 is not influenced by this value. The differ- 
ence of squared masses is of course very crucial. 
From k\ and k 2 one can determine k = k\ — k 2 . 
Then the oscillation frequency due to the inter- 
ference of the mass eigenstatcs m\ and m 2 is 
given by AE[ki, k, cos ip, m 2 ] as defined in Eq. (|7|). 
AE[k\, 0, cos (p, mi] is by definition zero. The an- 
gle ip can be put to zero, its absolute value should 
not be larger or very close to ir/2. As already 
mentioned this calculation has to be done with 
sufficient accuracy. At least 30 digits are recom- 
mended. Too large errors have catastrophic con- 
sequences. 

Detailed numbers for the EC-capture decay 
of 140 p r 58+ are shown in Table Q] For the sec- 
ond EC-capture reaction which was published 
in [2], the electron capture decay of the H-like 
ion 142 Pm 60+ to 142 Nd 60+ we get a prediction for 
Am\ x = 0.775(26) -lO^eV 2 . The corresponding 
numbers are listed in table [2] 

In conclusion, we are able to describe quanti- 
tatively the periodic modulations of the expected 
exponential decrease of the number of EC-decays 
per time which were observed in recent electron 
capture experiments in the Experimental Stor- 
age Ring (ESR) at GSI in Darmstadt [2] by neu- 
trino mixing. We give general formulae which re- 
late the modulation frequency with the difference 
Am 2 ; of squared masses of two neutrino mass 
eigenstates rrij and mi. With these formulae we 
get for both reactions where such modulations 
were observed, the nuclear K-shcll electron cap- 
tures of the H-like ions 140 p r 58+ and 142 Pm 60+ , 
values of Am^ which agree nicely with the global 
analysis of the solar-neutrino and KamLAND 
experimental data [6 . These general formulae 
which allow also to predict the modulation fre- 
quencies in further EC and bound beta decay ex- 
periments depend only on kinematical parameters 
of the process. We discuss that for the appearence 
of these modulations it is necessary to have exper- 
iments which produce a sufficient momentum un- 
certainty for the wave-functions (wave-packets) 
of the mother ions. The amplitude of the modu- 



li) 



lations is determined by the widths of the wave- 
packets, but the frequency depends on Am^ 
only. 

When this article was in the final stage of its 
preparation we realised that H. J. Lipkin has also 
published an article about the problem of neu- 
trino mixing in K-capture [TT]. In [T^] the same 
GSI-experiment [2] was explained by neutrino 
mixing using an ansatz for wave-packets in the 
final state only. 
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M d 

Q 


130319252000.000000000000000000000000000000000000000 
3388000.000000000000000000000000000000000000000 


TO 2 


0.100000000000000005551115123125782702118 
0.100380775051799640446407847251541936522 
0.000076299999999999997923709471603359588 


ki 
k 2 

K 


3387955.961051731594435668834899140341414454486 
3387955.961051731583175484621991824131728111711 
0.000000000011260184212907316209686342774 


Evi 
Ev 2 
Ev< 2 


3387955.961051733070251782586223414251214280672 
3387955.961051733070252075321238703166027210199 
3387955.961051733070252075321238703166027210199 


Edi 
E d i 

E 'd2 
E 'd2 ~ E d2 
E' m 2 

AE 


44.038948266929748217413776585748785719327 
44.038948266929747924678761296833972789800 
44.038948266929748510148791874663599621787 
0.000000000000000585470030577829626831986 
0.000000000000000000000000000000000486453 
0.000000000000000585470030577829626345533 



Tabic 1: Some quantities in eV, relevant for the explanation of the oscillation of the decay factor of the 
H-like ion 140 p r 58 + to 140 Cc 58+ . The calculations arc done for ip = according to Eq. Q and Eq. 0. 
toi is on-shell and to 2 off-shell. E V2 and Ed2 are the energies of m 2 and the daughter nucleus when the 
neutrino is emitted from the mother nucleus at rest. E' and E' d2 are the corresponding energies when 
TO2 is emitted from a mother nucleus moving with momentum —k against ki in order to get k<i = k\ — K. 



M d 

Q 


132321536600.000000000000000000000000000000000000000 
4870000.000000000000000000000000000000000000000 


TOl 
TO 2 

Am 2 21 


0.100000000000000005551115123125782702118 
0.100386752114011543561103471770365429653 
0.000077499999999999999986989573930173946 


h 
k 2 

K 


4869910.384855836089081365647620981219735553232 
4869910.384855836081124633213067424714490018121 
0.000000000007956732434553556505245535110 


Ev x 

EV2 

Ev' 2 


4869910.384855837115794303846383129612591300565 
4869910.384855837115794596682869979717867103615 
4869910.384855837115794596682869979717867103615 


Edi 

Ed2 
E 'd2 
E'd2 ~ Ed2 

E' m2 
AE 


89.615144162884205696153616870387408699434 
89.615144162884205403317130020282132896384 
89.615144162884205988990103720492684980937 
0.000000000000000585672973700210552084552 
0.000000000000000000000000000000000239217 
0.000000000000000585672973700210552084552 



Table 2: Same as Tabled] for the electron capture decay of the H-like ion 142 Pm 60+ to 142 Nd 1 
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Appendix 

Comments to Eq. ([2G]) 



1. Wc define the functions 

ss(x, t) ■ 

i(X,t) = 
I{t)-- 



sin[(x — Xo)t] sin[(a; + Xo)t] cos(2xot) — cos(2a;t) 



2(x 2 - xl) 



+x 



dxss(x, t) 



X 



dxss(x, t) 



(38) 
(39) 
(40) 



and get the analytical relation 



i(X,t) 



dxss(x, t) 



x 



dx 



sin[(x — xo)t] sin[(x + xo)t] 



= / dx 



cos(2a;o0 — cos(2a;i) 



2{x*-xl) 



2 r 

— < sin(2teo)r 
x I 

cos(2ta ) [ 



2t(X-x„) 



dx 
dx 



sini 

a; 
cosx 



2t(X+x ) 



dx- 



2t(X-x ) 



X 



2t(X+x ) 



X 

cos x X , 1 

da; 2atann — >. 

x - 1 



(41) 



x 



2. I(t) is periodic in t with to = 2xo and takes its maximal value at 2txo = 2mr + At such 
t-values the expression (|4Tj) simplifies to 



i(X,t) 



2t(X-x ) 



dx- 



XQ <- Jo 



2t(X+x ) 



dx- 



(42) 



The following diagram depicts i(X, t) at t = 8n + j for xq = 1. One can clearly see the 
influence of the two peaks of ss(x,t) at ±xq. 

1 . 75r 

1.5 




The values at these peaks ss(±Xo, t) oscillate with t and with an amplitude increasing propor- 
tional to t 



Ss(±Xq, t) = 



t sin(2a;oi) 
2xq 



(43) 



The widths of these peaks of ss(X,t) at X = ±xo and t — (n± t)— shrink proportional to 
1/t. 



13 



3. I(t) is oscillating with the same frequency like ss(±Xo,t). 

I(t) has a zero at t = n-^-. The following diagrams shows ss(X,t) for t = 8^-. 



1 


6 
4 
2 






-3 -2 


-2 
-4 
-6 


wwv\/u : 


2 



3 



J(i) has a maximum at i = (n + J) J 11 . The following diagram shows ss(X,t) for t = (8 + 1 ) ^- . 





10 - 
5 - 




-3 


« w v y y i 


||/VVWwww 


V| l|" — 2 


3 






-5 - 










-10 - 







I(t) has a zero at t = (n + \)^- The following diagram shows ss(X, t) for t = (8 + 5)^. 





1 10 

5 

IAMaaaa; 






-3 "- i 2 Wv ^- 


1 

-5 
-10 




IH^w^" """"" " 3 



has a minimum at f = (n+ §) The following diagram shows ss(X, t) for i = (8 + 





10 
5 

jjl/Wvwv 






-3 








y W v«— 2 - 


3 






-5 












-10 
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• Some modification of Eq. (|26|) 



sin[(fc - ki)L/2} sin[(fc - k 2 )L/2] |fci -fc 2 |z,/2s>2rr 
(k-h)/2 (k-k 2 )/2 Aj 

sin[(fc - fc 2 )L/2] sin[(fc - h)L/2] 



7T / dfc|<J(fc 



(k-k 2 )/2 



(k - ki)/2 



5{k-k 2 ))f{k) 



sin[(fci - k 2 )L/2] 
= n ( kl -k 2 )/2 VM + H 1 *)] 

Detailed derivation of Eq. ([27]) from Eq. (JUJ) 
1 d 



AbcW 



r d 3 fc 


f d 3 q 






(2nf/V 


L 



t 2 
dT(f\H W (T)\l MF ) 



Mq T d 



d 3 k 



d 3 q 



2 At J (2tt) 3 /V J (2tt) 3 /V 

d 3 K 



Mq T d 



y(47T 3 0-2) 3 / 4 

d 3 fc f d 3 q 



c d 3 Rc i{s -^- k)R J2 U ej dTCxp{iAEjt} 



2 dtj (2tt) 3 y (2tt) 3 
d 3 ^ 



Mq T d 



d 3 fc 



(47T 3 ct2)3/4 
d 3 K 



_K_ 2 3 

2CT 



II E ^ / d - exp{iAE jt } 



Kj —qj—k 
i=l 2 



2 dt J (27r) 3 7 (47T 3 Cr2)3/2 

x I y^Uej [ dTexp{iAE 3 t} 



dVe 



.5 

n[ 



dq . s in( K <-^- k > L) sm( K '- q *- ki 



(44) 



(45) 



D 



2ir K,- qi ~k t 

2 



k'- —qi—ki 



Mq T d 



d 3 fc 



2 dt./ (2tt) 3 7 (47r 3 a2)3/2 



d 3 A 



k 2 + k' 2 3 



dVe 



m 



sm( — 5—^) 



Mq T d 



d 3 k 



2 dt J (2tt) 3 7 (47r 3 a2)3/ 



d 3 K 



-5 

n[ 



I Vj7 ej / dTexp{iA^t} 
^C4j / dT exp{iA£^} 
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• Detailed derivation of Eq. (|29|) 



df 



' e 3 AE_ 
3 2 



. AEjt 

sin — ^- . a-Bj 
— e 2 



df 



sin — sm 

A~E~ 

2 



AS, -AE, 



2 



• 2 AEjt 

sm 2- 



dUd*^^ [^] 2 ■ E ••' -'"^i — a|— c ) 



d 
d7 



• AE,t ■ AEjt 

sm — ^ sm J 



^ ' AEi 
2 



Agj 
2 



S(AEi) 



(AEa - AS;)* 



} 



AE, 
2 



3 jjH 



■ AEit AEit ■ AEjt AEjt 

sm ±±f± cos sm — ^- cos — ^- 



AEj 
2 



5(A^)]} 



2tt^ L/£ <5(A£,) + 2wJ2 u ejUei cos(AEit)S(AEj) 

2nJ2HAE 1 )^ ] + J2 UejUei cos(A^t) 
i K¥>) 



(46) 



Reduction of the diagonal terms of Eq. (|27[) to Eq. (|33[) 



AecW 



Mq T d 
2 df 



Mq T d 



d 3 fc 
(2^)3 

d 3 fc 



d 3 , 



(47r 3 C72)3/2 
d 3 K 



m 

i=l 
3 



d«; 



,sin(^2^L) 



2 . ;2 
' 2a'i 



T, u > 



(2tt) 3 
d 3 fc 

2 dt J (2tt^ 



M 2 G t 



M GT 



M GT 



d 3 fc 



(2tt) 3 J (ira 2 K ) 3 / 2 
d 3 fc f d 3 K 



(^)3/2 

d 3 K -4 d X - 

c °% — > C7, 



dr exp{iA_Ej (fc, k, ip)t} 



AEit 

sm — j*- . ae 3 - 
2 e 2 



Agj 
2 



(2tt) 3 7 (7ra2)3/2 



M, 



GT 



d 3 fc 

(2tt) 3 



n 



2 

dr exp{iA£ , i3 -t} 



/ ^ / jj^sj* n [ / ^ 2 * - «o e " 2ct2 ' ] e ^ 1 1 
d 3 ^ 



,t 2 

dr exp{iAEjt} 



2tt5{AE) E C 7 , 



d3fc * rA/n _ M gt 4ttQ 2 _ Mq T Q 2 3 2 
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Reduction of a non-diagonal term of Eq. (f27|) to Eq. (|34l) 



d 3 fc -q [JdKi fdn'i 



K +K • / Ki— K, 



2 7 (2tt) 3 JLi L(47r 3 (j2)i/2 w< -< 

2 — 1 2 

(.Ue-jUei^- f dTexp{iAEj(k,K,',(p')t} [ dTcxp{—iAEit(k,n,tp)} 
dt J Jq 



f d 3 fc Yl\ l^I^j e -^f 



sin( K ' „ Ki L) 



2 J (27r)3llL(47r8^)V2 
x U e jU e i2-!r8(AEj(k,K,',ip'))coa(AEi(k,K,ip)t) = 



M GT 



d 3 fc 



n[ 



d« l 
"2T 



4crj 



■sin(f-L) 



x UejUefi-n-SiAEjik, 0, </?)) cos(A^(/c, k, </?)t) 
Folding of Gaussians 



d 3 fc 



d 3 g £ 2 +g 3 — 

e ~^~f(k-q) 



(27T(J 2 ) 3 / 2 7 (2trj 2 ) 3 / 2 
Integration of plane waves 

3 

J i=l 

s i n (9i±kL) 



d 3 /fc 



( 47r(7 2)3/2 



7(fe) 



d 3 <7 



a -i(g i + fc i )L/2 _ -i( 9l +fc i )L/2 



-i(qi + h) 



i=i 



q. + kj 
2 



JJ [27ri / d 9l c% + fe) = (2ttL) 3 

i=l 
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